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Abstract. Some standard examples of Hamiltonian systems that are integrable 
by classical means are cast within the framework of isospectral flows in loop alge- 
^ ! bras. These include: the Neumann oscillator, the cubically nonlinear Schrodinger 

^ . systems and the sine-Gordon equation. Each system has an associated invari- 

ant spectral curve and may be integrated via the Liouville-Arnold technique. 
^ ■ The linearizing map is the Abel map to the associated Jacobi variety, which is 

^ ■ deduced through separation of variables in hyperellipsoidal coordinates. More 

generally, a family of moment maps is derived, embedding certain finite dimen- 
sional symplectic manifolds, which arise through Hamiltonian reduction of sym- 
plectic vector spaces, into rational coadjoint orbits of loop algebras g"*" C g[(r)+. 
Integrable Hamiltonians are obtained by restriction of elements of the ring of 
spectral invariants to the image of these moment maps; the isospectral property 
^ i follows from the Adler-Kostant-Symes theorem. The structure of the generic 

5^ ! spectral curves arising through the moment map construction is examined. Spec- 

tral Darhoux coordinates are introduced on rational coadjoint orbits in 0[(r)+*, 
and these are shown to generalize the hyperellipsoidal coordinates encountered 
in the previous examples. Their relation to the usual algebro-geometric data, 
consisting of linear flows of line bundles over the spectral curves, is given. Apply- 
ing the Liouville-Arnold integration technique, the Liouville generating function 
is expressed in completely separated form as an abelian integral, implying the 
Abel map linearization in the general case. 
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1 Background Material and Examples 

In this first section, we shall examine several examples of integrable Hamiltonian 
systems that may be represented by isospectral flows on coadjoint orbits of 
loop algebras. In each case, the flow may be linearized through the classical 
Liouville-Arnold integration technique. An explicit linearization of the flows 
can be made in terms of abelian integrals associated to an invariant spectral 
curve associated to the system. A key element in the integration is the fact that 
a complete separation of variables occurs within a suitably defined coordinate 
system - essentially, hyperellipsoidal coordinates, or some generalization thereof. 
A general theory will then be developed in subsequent sections, based essentially 
on moment map embeddings of finite dimensional symplectic vector spaces, or 
Hamiltonian quotients thereof, into the dual space of certain loop algebras, the 
image consisting of orbits whose elements are rational in the loop parameter. 

The origins of this approach may be found in the works of Moser [Mo] on 
integrable systems on quadrics, Adler and van Moerbeke [AvM] on linearization 
of isospectral flows in loop algebras and the general algebro-geometric integration 
techniques of Dubrovin, Krichever and Novikov [KN], [Du]. The theory of 
moment map embeddings in loop algebras is developed in [AHP], [AHHl]. 
Its relation to algebro-geometric integration techniques is described in [AHH2] , 
and the use of "spectral Darboux coordinates" in the general Liouville-Arnold 
integration method in loop algebras is developed in [AHH3]. Some detailed 
examples and earlier overviews of this approach may be found in [AHH4], 
[AHH5]. The proofs of the theorems cited here may be found in [AHP], 
[AHH1-AHH3]. 

1.1. The Neumann Oscillator 

We begin with the Neumann oscillator system ([N], [Mo]), which consists of a 
point particle conflned to a sphere in R", subject to harmonic oscillator forces. 
The phase space is identiflable either with the cotangent bundle or the tangent 
bundle (the equivalence being via the metric): 

M = T*5"-^ = {(x, y) e M'^ X M'^ I x'^ x = 1, y^x = 0} C x M", (1.1) 

where x represents position and y momentum. The Hamiltonian is 



iy(x,y) = -[y^y + x^Ax], 



(1.2) 



where A is the diagonal nxn matrix 



A = diag(Q;i, . . . «„) e M 



nxn 



(1.3) 
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with distinct eigenvalues {ai}f=i,...n determining the oscillator constants. 
Equivalently, we may choose the Hamiltonian as: 

0(x,y) = i[(x^x)(y^y) +x^Ax- (x^y)2], (1.4) 

in view of the constraints 

x^x = 1, y^x = 0. (1.5) 
The unconstrained equations of motion for the Hamiltonian are: 

^ = (x^x)y - (x^y)x (1.6a) 

^ = -(y^y)x - Ax + (x^y)y. (1.6b) 

Since 

{0,x^x} = O, (1.7) 
it is convenient to interpret the relation 

n(x) := x^^x = 1, (1.8) 

alone as a first class constraint. The function n(x) generates the flow 

(x,y) ^ (x,y + tx), (1.9) 

and is invariant under the (p-Row. We may then apply Marsden-Weinstein re- 
duction, quotienting by the flow (1.9). The reduced manifold is identifled with 
a section of the orbits under (1.9) deflned by the other constraint: 

y^x = 0. (1.10) 

The integral curves for the constrained system are determined from those for 
the unconstrained system by othogonal projection: 

(x(t),y(t))/^ee^ (x(t),y(t)) constr. 

This lifts the projected flow on the Marsden-Weinstein reduced space n~^(l 
to one that is tangential to the section deflned by eq. (1.10). 
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Assuming the constants aj determining the oscillator strengths are distinct, 
the integration proceeds (cf. [Mo]) by introducing the Devaney-Uhlenbeck com- 
muting integrals 



which satisfy 



ai- aj 



n 



J^/j = x^x (1.13a) 

n 

Y^a,h = 2(t). (1.13b) 



Define the degree n — 1 polynomial P(A) by 

(114) 

o(A) ■ 4^A-ai' ^ ' 

where 

a(A) JJ(A - Qjj), 
i=l 

V{\) = Pn-lA"-^ + Pn-2A'^-2 + . . . + Po- (1.15) 

* cn— 1 



Then on T*S 

1 r 1 
P_i = --x x=--, 

Pn-2 = \<t>- (1.16) 

An equivalent set of commuting integrals consists of the coefficients of the poly- 
nomial (Pq) • • • 5 Pn-2}- The Liouville- Arnold tori T are the leaves of the La- 
grangian foliation defined by the level sets: 

Pi = Ci. (1.17) 

We now proceed to the linearization of the flows through the Liouville- 
Arnold method. First, introduce hyperellipsoidal coordinates {A^}^=i^...„_i and 
their conjugate momenta, {C^} i^i=i,...n-ij which are deflned by: 



n 



2 



£^A-a, " a(A) ^^'^^^^ 
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In terms of these, the canonical 1— form is: 

n n — 1 

e = ^yidxi\T*S"-^ ^^Ci^dXi^. (1.19) 

i=l /i=l 

Restricting this to T determines the differential of the Liouville generating func- 
tion 5': 



n— 1 n— 1 



1 V "'^^z^) 



which, upon integration, gives 



n-l „A 







This is seen to be an abelian integral on the (generically) genus g = n — 1 
hyperelliptic curve C defined by: 

+ a(X)ViX) = 0. (1.22) 
The linearizing coordinates conjugate to the invariants Pj are then: 

where, for = 2P„_2, 

6n-2 = 2, 6,=0, 3<n-2. (1.24) 

The map: 

(Ai,...A„_i)^(gi,...gn-i) (1.25) 

defined by eq. (1.23) is, up to normalization, the Abel map from the symmetric 
product S'^~^C to the Jacobi variety J{C) of C: 

A : S'^-^C ^ J{C) - C"-Vr, (1.26) 

where F = is the period lattice. 

We now turn to the interpretation of such systems as isospectral flows in a 
loop algebra (cf. [AHP], [AHH3], [AHH4]). Let 

J\f{X) = XY + J\fo{X), (1.27) 
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where 



2 I spn spn XiVi J - ' n no \ 

\ Z^i=l X-ai l^i=l X-ai / (1.28a) 

r:=(° 7). (1.28b) 

We define a map Ja IR"" x — > s[(2)+* to the dual space of the positive half of 
the loop algebra sl{2), relative to the standard splitting 5l(2) = sl(2)+* +s[(2)~* 
into holomorhic parts inside and outside a circle in the complex A-plane, 
containing the ctj's in its interior region: 

Ja: (x,y)^jVo(A)eJi:(2)+*. (1.29) 

This is a Poisson map with respect to the Lie poisson structure on 5l(2)+*. The 
Hamiltonian is then given by restriction of an elementary spectral invariant: 

0(x,y) = -tr(A^(A)2)o := / tr(A^(A)^)^, (1.30) 

and all the other invariants Pj may be similarly represented. The equations of 
motion are seen to be equivalent to the Lax equation: 

(1.31) 

where 

and {d4>{M)-^ signifies projection of the element d(f){J\f) e sl{2) to 5((2)+. This is 
an example of the Adler-Kostant-Symes (AKS) theorem (to be explained more 
fully in Section 2). The spectral invariants (elements of the AKS ring) are 
generated by the residues of the rational function 

,jm\^m^.^±^. (1.33) 

\ \ J a{\) A X-ai 

To see the relation with the standard algebro-geometric linearization meth- 
ods ([Du], [KN], [AHH2-AHH3]), we begin with the invariant spectral curve: 

det(£(A) - zh) = + a(A)P(A) = 0, (1.34) 

where 

C{X:=^N{\), z:=a{\)^, (1.35) 



•A/'o(A) ( "rr' ^r' 'r^ I ^ ^[(2)+* 
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and let 



A^(A,C):= (^-CIs), (1.36) 



where I2 is the 2x2 unit matrix. Then the columns of the matrix M.{X,Q of 
cofactors of A4{X,C) ^re the eigenvectors of A/'(A) {(X = eigenvalue) on C: 

(1 Y^"- XiVi _ > 1 v^n y'f -I- i \ 
2Z^^=lA-a. <^ 2 l^i=l X-a, ^ 2 \ ^^^.37) 

The hyperellipsoidal coordinates {Xn,(n}iJ.=i,...,n-i define the finite part of the 
zero- divisor: 

n-l 

V=Y,PiK,C^.)+P{ooi), (1.38) 

/x=l 

i.e., the zeros of a section of the bundle E ^ C dual to eigenvector line bundle. 
This bundle can be shown generically to have degree n, and thus to be an element 
of the Picard variety E G Pic"^. The Abel map then identifies the symmetric 
product S'^~^C with the Jacobi variety J{C) ~ Pic°. The linearity of the flow 
in Pic"^ follows from noting that the Lax equation 

^ = [dm)+,j^. (1-39) 

implies a linear exponential form for the transition function 

T(A,2,t) = exp((^^(A,^)t). (1.40) 



1.2 Nonlinear Schrodinger (NLS) Equation 

We now apply a similar analysis to the quasi-periodic solutions of the cubically 
nonlinear Schrodinger equation (cf. [PI], [AHP], [AHH4]) 

Uxx + \/-Twt = 2\u\^u. (1.41) 

Let 

/ • I g j p _ Zj \ 

m) ^ r ^ ^ ' (1-42) 

and let uj denote the standard symplectic form on C": 

n 

u> = i dz^ Adz = i dzj A dzj (1.43) 



We define the Poisson map: 



J : C" — > iu(l, 1)+* 
J:zi — >J\f{X), 



and let 



a(A) 



£(A) = ^A/'(A) = LoA"-i + LiA^-2 + • • • + 
A 



The spectral curve is defined by the characteristic equation 

det(£(A) - zl) =z^ + a{X)V{X) = 0, 

V{X) :^Po + PiX + --- + Pn-2A'^-^ 

and has genus g = n — 2. Choosing the AKS Hamiltonians: 



a(A) 



A tr(^(A)2) 



^A^ tr(A^(A)^) 



— -p2,n-3 



gives the Lax equations 



where 



^£(A) = [(di?.)+,£(A)] 
j^C{X) = [{dHt)+,C{X% 



{dH,)+ = ALo + Li 
{dHt)+ = A^Lo + ALi + L2. 



Choosing invariant constraints so that: 

" 2 
1/2 = i 







hip 






-\u\^ J 



u 

u or 



the compatibility conditions: 

d{dHx)+ d{dHt)+ 



dt 



dx 



+ [{dHx)+, {dHt)+\ = 
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reduce to the NLS equation (1.41). To obtain the quasi-periodic solutions, we 
introduce the spectral Darboux Coordinates {q, P, A^, Cix}^,=i,...n-2, analogous to 
the hyperellipsoidal coordinates above: 



n 



^ A - Q!i a(A) ^ ^ 



?:=ln(w), P:==.(Lo)22, (1.52c) 
Then the symplectic form may be expressed as 

n-2 

o; = ^ dA^ A dC/Lt + A dP. (1.53) 

As above, the spectral curve is invariant, and the coefficients of the charac- 
teristic polynomial generate a complete set of commuting integrals, so we may 
apply the Liouville- Arnold integration method. The coordinates (A^, C/i) defined 
by (1.52a,b) again give the finite part of the divisor of zeros of the eigenvectors of 
AA(A), while the remaining pair (g, P) are determined by the additional spectral 
data at A = oo. The Liouville generating function in this case becomes 



5= VCMciA^|p,=c.t.+9i'=y] / W ^dA + Plmx (1.54) 

and the linearizing coordinates conjugate to the Pj's are 

* = 11 = 5§i"7=|^ + ' = 0.-»-3. (1-55) 

which are abelian integrals of the first kind, and 

dS f^^ A--2dA Inw , 

^'^-' = W = -2}^X V-a(A)P2(A)"^ = ^-^" + ^-^^' ^'-''^ 

which is an abelian integral of the third kind, the integrand having simple poles 
at the two points (0x21,002) over \ — 00. For the Hamiltonian Hj; = —Pn-3, 
we have bi = -5i,n-3, Q = 0, while for Ht = -Pn-4, h = 0, q = -5i,n-3- An 
explicit formula for the function u{x, t) may be obtained in terms of the Riemann 
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theta function 9 associated to the spectral curve by applying the reciprocity 
theorem relating the two kinds of abelian integrals (cf. [AHH4]). 

u[x, t) - exp(gj - K exp(te + ct) ^(^(ooi,^) + tU + - K) ' ^^-^^^ 

where A : S^'-'^C i — > J{C) ~ C^-^/F is the Abel map, U, V G C^-^, &, c e C 
are obtained from the vectors with components Cj) on the RHS of eq. (1.55), 
(1.56) by applying the linear transformation that normalizes the abelian differ- 
entials in (1.55), and K is the Riemann constant. 



1.3 Sine-Gordon Equation 



As a last example, consider the sine-Gordon equation (cf. [HW] , [P2] , [A A] ) 



Let 



sin('u). (1.58) 



U{\):=XY + Uq{\), (1.59) 

where 

y=(j -;), A.„(A):=2A(_*(^ij 4^)), ,1.60) 
with 6(A), c(A) given by 



p 



i=l 



i=l 



and 7 e complex vectors with components {y'i, 7i}i=i,...p. (Here ctj, ctj, 
— Q!j, —oii are assumed to be distinct.) 

Define 

a(A):=n[(A^-«D(A' (1-62) 

The symplectic form on x 'CP is given by: 

p 

uj = 4^{d-fiAd(fii + d'yiAd(fii). (1.63) 

i=l 



Again, define a Poisson map: 

J: CP X CP — ^5u(2)+* 
J : {(p,^)^J\fo{X), 

where the twisted loop algebra: 

iu(2)+ c iu(2)+ c u(2)+ c gl{2, C)+ 

is defined as the fixed point set in sl(2, C)"*" under the involutions; 

ai : X(X) ^ X^(X) 

Let n — 2p, and 

£(A) := ^AfiX) 
A 

= a(A)F + LoA^"-! + LiA^"-^ + • • • + Lsn-i- 

The spectral curve is: 

det(£(A) - zl) = z^ + a{X)P{X) = 

P(A) = Po + A^Pi + • • • + A2^-2p„_i + A^". 

Choosing the AKS Hamiltonians: 

HdX) ^tr (^^(X(A) + AF)^)^ = -Po 
:= -itr (^(X(A) + AF)^)^ = P, 



n-l 



gives the Lax equations 



^£(A) = [A>C(A)] 



^£(A) = [P,/:(A)], 



where 



-A = dH^iAf)_ = j{L2n-i + a(O)y) 
B = dHr^{Af)+^Lo + XY. 
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Choosing the level set: 

Po=l^ (1.72) 

gives 

L2„_i + a(0)y = - \ j , (1.73) 

where 

= a(0)(c(0)-l), (1.74) 

with tt real. Then the compatibility conditions: 

dA dB , , , , , 

___ + [^,B]=0 (1.78) 

reduce to the Sine-Gordon equation 

Uxx - utt = sinu, (1-76) 

where 

^ = x + t, r] = x-t. (1.77) 

The quasi-periodic solutions are obtained as in the previous example. The 
unreduced spectral curve is defined by 

+a(A)P(A) = 0, (1.78) 

and is again hyperelliptic, with genus g — 2n — 1. Quotienting by the involution 

{z,X)^{z,-X) (1.79) 
gives a reduced curve with genus g = n — 1 defined by 

z'^ + a{E)P{E) =0, (1.80) 

where 

X'^:=E, a{Ey, = a{X), P(E) := P(A). (1.81) 
We also introduce the augmented curve, defined by 

z'^ + Ea{E)P{E) = 0, (1.82) 

of genus g = n, where 

z := zX. (1.83) 
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The spectral Darboiix coordinates are defined by 

c(E^) -1 = (1.84a) 

C^^ = 2b{E^), fi=l,...n, (1.84b) 

where 

b{E) := 6(A), c{E) := c(A). (1.85) 

These again are interpreted as zeros of the sections of the dual to the eigenvector 
line bundle associated to A/'(A). The symplectic form is then 

n 

u; = J2dE,,AdCi^ = -dd. (1.86) 
The Liouville generating function is 



r^^" / P(E) 

S{Po,...,Pn-i,Ei,...,En) = J2j^ \j-^^^dE, (1.87) 
giving rise to the Abel map linearization: 

^/.^i-^^o ^-Ea{E)P{E) 
= Ci + 25i,o^-2Si,n-iV, (l-88b) 

which only involves abelian integrals of the first kind on the augmented curve. 
In terms of theta functions, the solution may be expressed as 



u 



0(A(po, oo) - Ury - V4 - K) 

where U, V are again obtained from the coefficients on the RHS of (1.88) by 
applying the normalizing linear transformation to the abelian differentials ap- 
pearing in (1.88). Full details for this case may be found in [HW]. 

In the following two sections, a general approach to integrable systems 
is developed, yielding all the above results as particular cases, but allowing 
generalizations to more complex systems of higher rank. 
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2 Moment Map Embeddings in Loop Algebras 
2.1 Phase Space and Loop Group Action 

We begin by defining the generalized Moser space (cf. [AHP]) to be the sym- 
plectic vector space consisting of pairs (F, G) of rectangular N x r matrices: 

M = {{F, G) e M^''^ X M^''^} (2.1) 

with symplectic form: 

a; = tr dP'^ A dG. (2.2) 

The loop algebra, denoted q, consists of smooth maps from a circle S^, 
centred at the origin of the complex A-plane, into 0l(r), sl(r), or some subalgebra 
thereof. 

= qK"^) (or sl(r)) 
= {X(X) e Ql{r), XeS^ gCUoo}. (2.3) 

There is a natural splitting of g 

= 0++r, (2.4) 

as a vector space direct sum of the subalgebra 0+, consisting of elements X{X) 
admitting a holomorphic extension to the interior of S^, and 0~, consisting of 
elements admitting a holomorphic extension to the exterior, with normalization 
X{oo) = 0. We identify g as a dense subspace of its dual space g* through the 
pairing 

< X > := -i. / tr (/.(A)X(A)) ^, (2.5) 

2m J SI A 

e 0_ , X eQ+. 

Under this pairing, we have the identification 

(5+)*^5_, (r)*-5+, (2.6) 

where 

r = 0++0- (2.7) 

similarly represents a decomposition of into subspaces consisting of elements 
holomorphic inside and outside S'^, but with the normalization such that ele- 
ments e satisfy //(O) = (and hence the constant loops are included on 

0-)- 
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The loop group 6"'' is similarly defined to consist of smooth maps g : ^ 
Gl{r) which admit holomorphic extensions to the interior of S^: We define a 
Hamiltonian action: 

6+ : M — >M 
giX):iF,G)^iF„Gg), (2.8) 

where {Fg,Gg) are determined by the decomposition 

{A - XI)-^Fg-'{X) = {A- Xir^Fg + Fhoi (2.9a) 
g{X)G^{A - A/)-i = G^{A - XI)-' + Ghoi- (2.9b) 

Here A G M^'^ is some fixed N x N matrix, with eigenvalues in the interior 
of S', and (-Fhoi, Ghoi) denote the parts of the expressions on the left that are 
holomorphic in the interior of S'. This Hamiltonian action is generated by the 
equi variant moment map: 

: M ^ 5+ ~ 0_ 
J^{F, G) = XG'^iA - XI)-'F, (2.10) 

which is thus a Poisson map with respect to the Lie Poisson structure on 0+* . 
This map is not injective, its fibres being (generically) the orbits of the sub- 
group Ga '■= Stab(^) C Gl{N) acting by conjugation on A, and by the natural 
symplectic action on M. 

g : (F, G) {gF, {g^)-'), g e Ga C Gl{N) (2.11) 

The relevant phase space is therefore the quotient 

M/Ga^Qa, (2.12) 

which is identified with a finite dimensional Poisson subspace C g"*"* con- 
sisting of elements that are rational in the loop parameter A, with poles at the 
eigenvalues of A. 

2.2 Simplest Case 

We now consider the simplest case, where A is a diagonal matrix: 



A = diag(ai, ... ai, ... a^, ... a^, «„), 



(2.13) 
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possibly with multiple eigenvalues {(Xi}i=i^,,,n of multiplicity ki < r, all in the 
interior of S^. The matrices (F, G) are decomposed accordingly as 



F = 



\FnJ 



G = 



Gi 



\GnJ 



(2.14) 



where (Fj, Gi) are the kiXr dimensional blocks corresponding to the eigenvalues 
CKj. For this case, Nq{\) has only simple poles: 



GjF, 



Mo{\) = J^{F,G) = -\Y,T± 



i=l 



A — a,- 



A — a,- 



(2.15) 



with residue matrices Ni generically of rank: 

rk(Fi) = rk(G'i) = h 



The Ad* -diction for this case becomes: 



(2.16) 



(2.17) 



which can be identified with the ^(i*-action of the direct product group Gl{r) x 
• • • X Gl{r) (n times) on [^[(r)*]". 

The fibres of the map coincide with the orbits of the block diagonal 
subgroup: 



Ga = Stah{A) = Gl{ki) X Gl{k2) x • • • x Gl{kn) C Gl{N), 



under the action: 



(/ii, . . . /i„ . . . hn) : (F„ G,) ^ {h,Fi, {hJy^Gi) 



(2.18) 



(2.19) 



This is also a Hamiltonian action, generated by the "dual" moment map: 

Jh{F, G) := (FiGf , ....... FnGl) e (0l(A;i) x • • • x g^A;,))*. (2.20) 

The Ad*^^ orbits are then the level sets of the Casimir invariants: 



tr(F,Gf)^ k = l,...n l = l,...ki 



(2.21) 



17 



More generally, the image of the moment map is a Poisson submanifold 
of C 'q'^* consisting of elements of the form 

n Pi 

^oW = ^T.j:aZ^' (2-22) 

where pi is the dimension of the largest Jordan block of A corresponding to 
eigenvalue ctj. 

2.3 Dynamics: Isopectral AKS Flows 

The Hamiltonian flows to be considered are those generated by elements 
of the ring of Ad*-invariant functions restricted to the translate \Y -\-'qa 

of the subspace 'qa by a fixed element \Y G 0~*, where Y e 0l(r). (The latter 
is an infinitesimal character for 5~, since it annihilates the commutator of any 
pair of elements.) We denote the ring of elements so obtained by 

^lKS--=mi\xY+~,.^ (2-23) 
and refer to it as the AKS (Adler-Kostant-Symes) ring. 
Let 

M{\) = \Y ^Mq{\) eXY + QA- (2.24) 

Wc then have the fundamental theorem that underlies the integrability of the 
resulting Hamiltonian systems, the Adler-Kostant-Symes theorem: 

Theorem 2.1 (AKS). 

1. If H & ^AKS' Hamilton's equations are: 

Xh{N) = ^ = [{dH)^,^^\ = -[{dH)_,J\f] (2.25a) 

2. ifH^.H^ell^^, 

{Hi,H2} = 0. (2.25b) 

Thus, all the AKS flows commute, and arc generated by isospectral defor- 
mations determined by Lax equations of the form (2.25a). In fact, it may be 
shown ([RS], [AHP], [AHH2]) that on generic coadjoint orbits of the form 
(2.22), these systems are completely integrable; i.e., the elements of the Poisson 
commutative ring X^j^g generate a Lagrangian foliation. Since the map (2.10) 
with image consisting of elements of the form (2.22) is a Poisson map, and passes 
to the quotient Poisson space M/Ga to define an injective Poisson map, the same 
results may be applied to the puUback J"^ o H oi any Hamiltonian in the AKS 
ring JXks- 
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Corollary 2.2. The results of Theorem 2.1 remain valid if the Hamiltonians 
Hi, H2 are replaced by o Hi, o H2 on the space XY + 'qa, identified with 
M/Ga through 

= Ay + jVo (A) = Ay + Xi {F, G) . (2.26) 



2.4 Reductions 

To obtain interesting examples, one usually must reduce the generic systems 
described above in a manner that is consistent with the structure of the dy- 
namical equations. This generally consists of Hamiltonian symmetry reductions 
involving either continuous or discrete symmetry groups. (It may also involve 
symplectic, or more generally, Poisson constraints.) We briefly summarize the 
procedure for both types of symmetry reductions below. The discrete Hamilton- 
ian reduction procedure is described in greater detail in [HHM] ; the continuous, 
Marsden-Weinstein reduction is fairly standard [AM] . 



2.4.1 Discrete Reduction: 

We consider discrete groups generated by elements a either of finite order 
or generating compact orbits, which act on the space M by symplectic diffeo- 
morphisms, and as automorphisms of the loop algebra g""". Let 



a : M 

be such a symplectomorphism, and 



M 



the corresponding automorphism of g"^ , with dual Poisson map 



+* 



+* 



(2.27) 



(2.28a) 



(2.28b) 



We assume that the moment map J"* intertwines these two actions, so that the 
following diagram commutes 



M 



M 



J' 
+* 



(2.29) 



It follows that J"^ may be restricted to the fixed point sets 



CM, l+:=0+C0+, 



(2.30) 
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and its restriction defines a moment map from the fixed point set M^j to the 
dual space 

I+* := 0+* C 0+* (2.31a) 

of the subalgebra 

1+ := 5+ C 5+ (2.31b) 

of fixed elements under a^. The results of Theorem 2.1 and Corollary 2.2 may- 
then be applied on the reduced spaces, provided the Hamiltonians in the ring 
IXks chosen to be invariant under the symmetry a*. 

2.4.2 Continuous Hamiltonian Reduction 

All the Hamiltonians in the ring IXks invariant under the Hamiltonian 
group action given by conjugation of AfiX) by A-independent elements in the 
stability subgroup of Y: 

Gy := Stab(F) C Gl{r), Qy := stab(F) C Qi{r). (2.32) 

This action is generated by a moment map Jy, given by the leading term A^o of 
A^o, restricted to Qy 

Jy:=No\s^, (2.33) 

where 

J^=A'o(A) = Ao + AiA-i + .... (2.34) 

Since the elements of ^Xks invariant, Jy is conserved under all the AKS 

fiows. Fixing a level set 

JY = fio^ 0y, (2.35) 

which we assume to be a regular value of Jy, and restricting to the coadjoint 
orbit CaToCA) C (0"*")*, the reduced space is 

0,ed■.= JY\^^o)/Go, (2.36) 

where Gq C Gy denotes the stability subgroup of hq. The reduced Hamiltonians 
i^red on Cred Sixe then given by 

HredOT^ = H\j-i(^^^-^, HeIXksi (2.37) 

where 

TT : Jy\fio) Jy'M/Go (2.38) 
denotes the projection map. 
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2.5 Examples 



We now indicate how the loop algebra formulation of examples like those of 
Section 1 is obtained from the general scheme described above. 



2.5.1 Neumann Oscillator (and similar examples) in 5l(2, R) 



Consider the case r = 2,ki = l,n = N. A discrete antilinear involution 
gives the reality conditions: 



CKj = CKi, F — F, G = G 



(2.39) 



reducing 0[(2, C)* to 0[(2, R)*. The stabilizer oi A = diag (q!i, . . . , an) consists 
of the diagonal subgroup Ga = {diag((ii, . . . , dn) C Gl{n)} acting as 



Ga--M — > M 



{di, ...dn) : 



\ ; / 



d^F, 

V ; / 



d-^G. 



\ 



(2.40) 



where Fj and Gj are just 2-component row vectors. The moment map generating 
this action is just 



JAiF,G) = {F^G^,...,FnGl)e 



(2.41) 



which coincides with the traces of the 2x2 residue matrices Ni in (2.15). Choos- 
ing the zero level set for these, Marsden-Weinstein reduction is equivalent to the 
subgroup reduction 0t(2)+ D sl(2)+. Choosing an appropriate symplectic section 
gives the reduced parametrization: 



-L(x,y), G=i=(y,-x) 
x,y eM". 



(2.42) 



The reduced symplectic form becomes 

(jj = dx^ A dy. 



(2.43) 



The reduced moment map is 



J\fo{X) = J''{F,G) 



A 



En Xj-iji 
i=l A— Of; 

^1=1 X—ai 



En xj-y. 
i=l 



(2.44) 



A— a; 
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Viewing this as defined on the symplectic vector space R"' x R", there 
is a residual fibration generated by the finite group (Z2)" of refiections in the 
coordinate hyperplanes. The sl{2)~ character XY may be expressed as: 

XY = X^l _^Je5l(2)-*, (2.45) 

and the resulting AKS flows involve isospectral deformations of elements of the 
form: 

-a)+-^o{X). (2.46) 

The Hamiltonians are chosen, as usual, from the AKS ring X^j^g(sl(2)*). 
In addition to the symmetry reductions already implemented, it is possible to 
impose further symplectic constraints of the form 

/(x, y) = 0, ^(x, y) = 0, {/, g} ^ 0, (2.47) 

and apply the standard methods for constrained systems. (Provided one of 
these functions is in the Poisson commutative ring Xj^p^g(s[(2)*), the constrained 
Hamiltonians will still commute.) The particular case of the above with 

= 0, b = -^, c = 0, (2.48a) 
/:=x^x-l = 0, ^:=y^x=0, (2.48b) 

and Hamiltonian (1.30), gives the Neumann oscillator system. The invariant 
spectral curve is of the form 

det(£(A) - ^I) = ^2 + a{\)V{\) = 0, (2.49) 

where 

C{X) := ^A^(A), (2.50) 

and 'P(A) is generally a polynomial of degree n — 1 or n, depending on whether 
+ be vanishes or not. 

2.5.2 The NLS Equation: Reduction to su(l, 1)+ 

Again, wc choose r = 2, ki = l,n — N. Similarly to the previous example, 
the zero moment map reduction under 

Ga = Stab{A) = X X • • • X (2.51) 
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is equivalent the subgroup reduction 0l(2, C)+ D sl(2, €)+. Choosing a suitable 
symplectic section gives the parametrization 



F=-L(.,w). G=^(w,-.), 



z, w e C". 



The further reality conditions 



(2.52) 



-^(z,iz), G=-^(-iz,z) 



(2.53) 



give the discrete reduction sl(2, €)+* D su(l, 1)+* as the fixed point set under 
an antilinear involution. On this real subspace, the symplectic form becomes 



uj — idz^ A dz 



(2.54) 



and the reduced moment map has the form. 



M)(A) = jV,g) = M '^^=^ 

2 I _ v^J^ 

\ Z^j=l X-a 



J = l X-aj 

2 




(2.55) 



In this case, we choose the character XY to vanish, so M{X) coincides with A/o(A). 
The commuting flows are generated by the pair of commuting Hamiltonians: 



Hx = 



a(A) 



A tr(AA(A)2) 



^A^ tr(AA(A)^)' 



e2^Ks(5u(l,l)+*) 
e2^Ks(5u(l,l)+*). 



(2.56a) 
(2.56b) 



Further invariant constraints are added to ensure that the leading terms of the 
polynomial matrix (1.45) have the form given in eq. (1.49). Defining C{X) again 
as in eq. (2.50), the resulting invariant spectral curve again has the form 



det(£(A) - zl) =z'^ + a(A)P(A) = 0, 
£(A) := ^M{X), 



(2.57a) 
(2.57b) 



where V{X) now is of degree n — 2. The Lax form (1.48a,b) of Hamilton's 
equations then follows from the AKS theorem, and the compatibility condition 
(1.51) gives the NLS equation (1.41). 
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2.5.3 Higher Rank Case. Two Component Coupled NLS System: 
Reduction to 511(1,2)+ 

As an illustration of a system described by an algebra of higher rank, we 
consider the case of the coupled two component cubically nonlinear Schrodinger 
equation {viz. [AHP], [AHH2], [AHH3]): 

mt + tta;a; = 2tt(|ttp + |vp) (2.58a) 
ivt + Va;x = MH'^ + \v\^)- (2.58b) 



In this case, we take r = 3 and ki = 1 for all z, so n = N. The process 
of discrete and continuous symmetry reduction is applied analogously to the 
preceding case, giving the sequence £|[(3, C)+* D s((3,C)+* D iu(l,2)+*. The 
reduced form of the resulting pair of n x 3 matrices (F, G) is 



F={p,ri,C), G={p,-ri,-0, 



(2.59) 



where rj, ( E C" is a pair of complex n-component column vectors and p e 
is a real column vector with components 



Pi = + ICil^, i = 1, 

The reduced symplectic form is 



. n. 



u = i{drj'^ Adr] + dC A dC), 



(2.60) 



(2.61) 



so the components (ry^, Ci)i=i,...n and their complex conjugates provide a canon- 
ical coordinate system on Oj^f^. The reduced moment map has the form 



M„(X) = J-\F,G) 



A - a,; 




ViPi 

-dm 



dp 
-IC 



(2.62) 



so the coadjoint orbit may be identified with x C". Again we choose the 
character XY to vanish, so J\f{\) coincides with A/o(A). As for the single compo- 
nent NLS equation, the commuting pair of Hamiltonians for the two component 
CNLS case is chosen to be 



a{\) 



eT(iu(l,2)+*) 



^A^ tr(A^(A)^) 



eX(iu(l,2)+*) 



(2.63a) 
(2.63b) 
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Defining, as before, 



jC(X) := ^XoiX) = LoA-i + LiA"-2 + • • • + L„_i, (2.64) 

further invariant constraints must also be imposed, implying that the leading 
terms are of the form (cf. [AHP]): 








u 






■ 









V V 





0/ 




(2.65) 



The Lax equations generated by the Hamiltonians (2.63a,b) have the same 
form as eqs. (1.48a,b), and the compatibility conditions (1-51) are equivalent to 
the CNLS system (2.58a,b). The invariant spectral curve in this case is a three 
sheeted branched covering of determined by an equation of the form 

det(£(A) - zl) = z^ + a{\)zV{\) + a{\fQ{\) = 0, (2.66) 

where V{X) and Q(A) are polynomials of degrees n — 2 and n — 3, respectively. 



3. Spectral Darboux Coordinates and Liouville- 
Arnold Integration 

In this section, the general method of linearization of AKS flows in rational 
coadjoint orbits will be explained. For simplicity, the spectral properties of the 
matrix A will be chosen as in Section 2.2, but the method is equally valid in 
the more general case (see [AHH3]). It consists of two steps. First, a suitable 
generalization of the hyperellipsoidal coordinates encountered in the examples 
of Section 1 is introduced, the spectral Darboux coordinates (Theorem 3.2) asso- 
ciated to the invariant spectral curve C. These consist of families of canonical 
coordinates on coadjoint orbits Oj^^ of the type discussed in the preceding sec- 
tion, which are naturally associated to the spectral data of the matrix A/'(A). 
The second step consists of using a Liouville generating function S to compute 
the canonical transformation to coordinates conjugate to the spectral invariants, 
in which the flow becomes linear. It turns out that for all Hamiltonians in the 
AKS ring IXks ^^^^ generating function, deflned with respect to the natural 
isospectral Lagrangian foliation of Oj^/^, may be expressed within the spectral 
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Darboiix coordinate system in completely separated form (Theorem 3.3). It fol- 
lows from the construction that this transformation is given in terms of abelian 
integrals, showing that, in the general case, the Abel map yields a linearization 
of the flows on the Jacobi variety J{C) of the spectral curve. One thus arrives at 
the algebro-geometric linearization results [viz. [Du], [KN], [AvM]) entirely 
through classical Hamiltonian methods. 

3.1 Phase Space and Group Actions 

In the following, the phase space will initially be thought of as a coadjoint orbit 
CaTo within the image of a moment map of the type introduced in Section 2. 

J^:M — (3.1a) 
: (F, G) I — > \G^{A - WY^F. (3.1b) 

The image defines a finite dimensional Poisson submanifold 

Im( J^) := C 0+* (3.2) 

which, in the simplest case, consists of elements of the form 

0^ = {M)(A) = Aj]— ^}, (3.3) 



where the ranks {/ci}i=i,...n of the residue matrices Ni coincide with the multi- 
plicities of the eigenvalues {Q;i}i=i,...n of the diagonal N x N matrix 

A = diai.g{ai,...,ai,...an). (3.4) 

The coadjoint action of the loop group 0'^ on Qa in this case becomes equiv- 
alent to the coadjoint action of {Gl{r))'^ on (^[(r)*)'^, obtained by evaluating the 
group element ^(A) e <3'^ at the poles X = af. 

g : {Ni} ^ {g{ai)Nig{ai)-'}. (3.5) 

It follows that the Ad~ -orbits are determined as simultaneous level sets of the 
Casimir invariants of the separate residue matrices Ni under this action: 

" N- 

OMo = {A V -^-^ I triVi = au 1 = 1,... h). (3.6) 
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The equations of motion induced by any element of the AKS ring e ^aks 
have the Lax form: 

^ = lii.(MW,M], (3.7) 

where 

^=*lAr+^. (3-8) 

is the restriction of the ^d* -invariant element $ e to the subspace con- 

sisting of elements of the form 

Af{X) = XY + Afo{X), AfoeQA. (3.9) 
Define the 5[(r)-valued polynomial 

£{X) := ^AA(A) = a{X)Y + LoA^-^ + • • • + L^-i, (3.10) 
A 

where 

n 

a{X):=Y[{X-ai) (3.11) 

i=l 

is the minimal polynomial of A. This satisfies the equivalent Lax equation 

^^mAfU,C]. (3.12) 

The invariant spectral curve Cq is then determined by the characteristic equation 

det (£(A) - zir) = (3.13) 

which, after suitable compactification, is viewed as an r-fold branched cover of 
P^, possibly having singularities over the points A = due to the r — ki fold 
multiplicity of zero eigenvalues. Other singularities could, of course, also occur, 
but for simplicity we again place ourselves in a "generic" situation in order to be 
able to give the main results in as explicit form as possible, and therefore exclude 
this possibility. The essential results remain valid without such simplifying as- 
sumptions, but explicit formulae for the spectral polynomial, genus, dimensions 
of orbits, and form of the abelian difi'erentials must be modified accordingly. 

We assume henceforth, for simplicity, that the spectral curves Cq have no 
singularities other than those that arise over {A = ai}i=i,...n5 if /ci < r — 1, 
due to the multiple zero eigenvalues of the residue matrices A^^. This implies 
in particular that the A^i's, while not necessarily diagonalizable, must lie on 
orbits that have the same dimensions as the diagonalizable orbits whose nonzero 
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eigenvalues are distinct; namely, ki{2r — ki — l). We also assume that one of two 
conditions holds, which excludes further singularities over A = oo: 

Case (a) Y = and Lq lies on a Gl{r) orbit of the same dimension (r(r — 1)) 
as those with simple spectrum. 

Case (b) Y and lies on a Gl{r) orbit of the same dimension as those with 
simple spectrum. 

Remark. An effect of this assumption is to eliminate from consideration the 
example 2.5.3, which has singularities over A = oo. However, this case may also 
be dealt with {viz. [AHH3]), by imposing a further set of symplectic constraints 
defining a "generic" deformation class of admissible spectral curves. 

In some cases, it is not the orbit Oj^^ itself that is the relevant phase space, 
but its reduction under the Hamiltonian action consisting of conjugation by the 
stability subgroup Gy = Stab (Y) C Gl{r): 

g : Mo{\) ^ gMo{\)g-\ geGy. (3.14) 
The corresponding moment map is just the leading term in Ao(A): 

n 

J(AAo):=Lo = 5^iVi, (3.15) 

restricted to the subalgebra Qy '■= stab(F) C 0^(r). Another case of interest, 
particularly when Y = 0, consists of restricting to a symplectic submanifold 
Ofj-^^ C 0_\f^ determined by the zero level set of the components of Lq within the 
annihilator of a Cartan subalgebra. (This is symplectic at regular elements Lq.) 
For future reference, we list the various subcases of interest. 

Case (a) Y = 0. In this case, all the elements e AKS ring are 

invariant under the full Gl{r) action (3.14), and all components of Lq are con- 
served. We may therefore reduce by the full group Gl{r) or any of its subgroups. 
The two subcases of greatest interest are: 

Case (a.l) Complete reduction at a regular point Lq = //q G sK'^)*- The reduced 
manifold is then 

= J-'M/Gq, (3.16) 
where Gq C Gl{r) is the stabilizer of hq. The dimension of the reduced orbit is: 

dim O'j^^ = dim O^r, - (r - l)(r + 2). (3.17) 

Case (a.2) Symplectic invariant manifold. We take the zero level set of all com- 
ponents of Lq in the annihilator of a Cartan subalgebra, (e.g., we choose Lq to 
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be diagonal). Denote this submanifold, which is symplectic at all regular values 
of Lq, as: 

Of;^ := {Nq e Ou^\Lq e r (Cartan subalgebra)}. (3.18) 
Its dimension is 

dim Of;^ = dim 0^f, - r{r - 1) = dim O'j^^ + 2(r - 1). (3.19) 

Case (b) Y 0. In this case, the elements e -^\'^+*+\y AKS ring are 

only invariant under the action of the stabilizer Gy = Stab(y) C Gl{r). Two 
cases of special interest arise: 

Case (b.l) The full orbit Oj\f^ (i.e., no reduction). 

Case (b.2) The reduction of Cjv'o under the full stabilizer Gy C Gl{r) of a 
regular element Y e gl{r)^ taken at a value I/olgy = A*o £ 0r- The group Gy is 
a maximal abelian subgroup with r — 1 dimensional orbits and Gq = Gy. The 
reduced orbit is denoted 

O^^r'^Jl^iM/Gy, (3.20) 

and has dimension 

dim cX^f = dim CaAo - 2(r - 1). (3.21) 



3.2 Structure of the Spectral Curve 

The affine part of the spectral curve Co is determined by the characteristic equa- 
tion (3.13). Taking into account the ranks of the residue matrices {A^i}i=i,...n 
in (3.3), we see that the characteristic polynomial has the general form 

V{X, z) = det (£(A) - zir) 

r 

= i-zY + z^-'V,iX) + (3.22) 

where 

n 

A, (A) := JJ(A - a^)'^ax(Oj-fcO^ rank £(«,) = ki. (3.23) 

This shows that near A ~ oo, we have 

_ ( n if y = (case (a)) , , 

zr^X"^, m:=l 3.24 
\n-l if r ^0 (case (b)). ^ ^ 
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This suggests a compactification, not within P^, but rather in the total space of 
a line bundle over ¥^ = Uo[J Uoo (where Uq^ Uoo denote the open disks obtained 
by deleting A = oo and A = 0, respectively), with coordinate pairs (A, z) over Uq 
and (A, 'z) over Uoo related by 

{\z)^(X=j,z=^) (over C/onC/oo). (3.25) 

This is just the total space T of the bundle 0{m) whose sheaf of sec- 

tions consists of homogeneous functions of degree m. The transformation (3.25) 
extends the affine curve Co defined by (3.13) over A = oo, defining the compact- 
ification: 

Co^C^T. (3.26) 

The possible spectral curves so arising are branched r-sheeted covers of 
P^, which within any given orbit of type (3.3), have ^-values over each A = 
that are fixed (being Casimir invariants of the coadjoint action (3.5)). Of these, 
there are ki nonsingular points (A — ai,z — Cia)a=i,...ki corresponding to the 
nonzero eigenvalues of C{ai), and the point (A = ai,z = 0), which generically 
is an r — ki-fold ordinary singular point corresponding to the r — /c^-fold zero 
eigenvalue. Figure 3.1 gives a visualization of the spectral curves C, embedded 
in T, as branched coverings of P-*^, constrained to pass through these points. 



Figure 3.1 

The detailed structure may be expressed more precisely by writing the 
form of the characteristic polynomial V{X,z) for any spectral curve C in a neigh- 
borhood of a given curve Cr as a perturbation of the characteristic polynomial 
Vr{\z) defining Cr. 
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Proposition 3.1 ([AHH3]). In a neighborhood of the point Mr G Oj^^ with 
characteristic polynomial Vr{X, z), the characteristic polynomial has the form: 

r 

V{X,z) = VR{X,z)+aiX)J2a,{X)p,{X)z^-^, (3.27) 

where 



Oj.(A) = Y[{X - (o.J-'^i-i) (3.28) 

i=l 

p,(A)=:^P,aA«, (3.29) 



a=0 



and {Pj{X)}j=2,...r ore polynomials of degree : 

, . \ di-j if Y = 
{ dj if Y ^0 

n 

dj = min{j — l,ki). (3.30b) 

The number of independent spectral parameters {Pja}, (a — 0, . . . Sj + n — m — 
1, j — 2, . . .r) is thus: 

d = g + r-l, (3.31) 



where 



1 1 " 

^= -(r - l)(mr - 2) - - ^(r - h){r - h - 1) (3.32) 

i=l 

is the genus of the (partially) desingularized spectral curve C obtained by sepa- 
rating branches at {ai, 0}. In a neighborhood of any generic point on Oj\f^, these 
spectral invariants are all independent. 

The complete integrability of the systems under consideration on the var- 
ious coadjoint orbits Oj\f^, and the reductions C^^, O^^J"^ and symplectic sub- 
manifolds Off^ thereof, may be seen from the following Table of generic dimen- 
sions for the various cases discussed above. (Note that the P^a's referred to do 
not include the leading terms of the polynomials Pj{X) in eq. (3.29).) Recall that 
the value of the genus g of the desingularized curve C given in Proposition 3.1 
depends on the value m, which is different in the cases (a) and (b): 

_ _ r n — 1 for case (a) 

g = g{m), m = \ (3.33) 
n tor case (b). 
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Table of Dimensions 

Case Dimension i^PiaS #-Pi's 

9 r-1 

£ 

9 r-1 
9 

Here, the notation {Pi}i=2,...r is used to denote the components of Lq 
evaluated on a basis of the relevent Cartan subalgebra (not including the trivial 
central element, which is a Casimir invariant, and hence constant on orbits). 
These are also elements of the ring X^Kg, corresponding to the leading terms in 
the polynomials Pj{X) in eq. (3.29) for case (a) , and the next to leading terms for 
case (b) (the leading terms being constant in the latter case), but they are listed 
separately since, in cases (a.l) and (b.2), they are fixed through the Marsden- 
Wcinstcin reduction procedure, and hence do not contribute to the number of 
independent invariants on the reduced spaces. Moreover, these elements enter 
again in Section 3.4 when defining the spectral Darboux coordinates for these two 
cases. It follows from the dimensions in the Table and the independence of the 
commuting invariants that cases (a.l), (a.2), (b.l) and (b.2) all give completely 
integrable systems. 

3.3 Spectral Lagrangian Foliation 

The Lagrangian foliation given by fixing simultaneous level sets of the invariants 
in the ring ^Xks spectral invariants (i.e., fixing the spectral curve C) is depicted 
below in Figure 3.2 for the various cases discussed above. 



Of,^ or O^, (O^J or O^'") 



Liouvillc- Arnold 
Tori T — > 
(isospectral leaves) 
dim T = 
g + r-1 (org) 











Admissable curves C G T 
dim C = g + r — 1 (or g) 

Figure 3.2 



We may summarize the relevant spectral data associated to each element 
J\f e Ay + qa as follows: 



(a) dim O^ro = 2g + {r - l){r + 2) 
(a.l) dim O^j^^ = 2g 

(a.2) dimO^^ =2(5^ + r-l) 

(b.l) dimOj^,=2{g + r-l) 

(b.2) dim O^^'"' = 2g 
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• A spectral curve C (r-fold branched cover of P-*^) defined by the characteristic 
equation: 

V{X, z) = det(£(A) - zlr) = 0, (3.34) 

(after suitable compactification and desingularization) . The r — 1 points over A = 
oo are determined by the leading terms of the polynomials pj (A) of Proposition 
3.1. 

• An eigenvector subspace: [V{X, z)] C at each point in C which, by our 
genericity assumptions, is one dimensional. 

Together, these determine an eigenvector line bundle E ^ C, and its cor- 
responding dual bundle E C. The latter may be shown {viz. [AHH2], 
[AHH3]) to be generally of degree 

deg{E)=g + r-l, (3.35) 

and hence, by the Riemann-Roch theorem, to have an r-dimensional space of 
sections. Conversely, it turns out that this data is sufficient to reconstruct the 
matrix £(A) (and hence Af{X)) up to conjugation by an element of Gl{r); that 

y red. 

is, it is sufficient to determine the projected point in the reduced orbit Oj^^ or 
C^J, but not the element M{X) itself. 

In order to reconstruct the element A/'(A), it is necessary to add some 
further spectral data, consisting of a framing at A = oo; that is, a basis of 
sections {ai G H^{C, E)}i=i^,,,r of the bundle E ^ C, chosen to vanish, e.g. 
at all but one of the r points {ooi}i=i^,,,r over A = oo (for the case where the 
Cartan subalgebra in question consists of the diagonal matrices). 

c7j(oOj)=0, iy^j. (3.36) 

This adds r — 1 dimensions to the fibres, (since framings related by {ai = nai} 
are equivalent). Furthermore, the spectrum over A = oo in the class of admissible 
spectral curves must be left undetermined, adding r — 1 dimensions to the base 
space in Figure 3.2. 

More generally, it is insufficient to just consider line bundles, since this ex- 
cludes the possibility of degeneracy in the spectrum and further singular points. 
The appropriate generalization consists of a coherent sheaf defined by the exact 
sequence 

0{-m)®' o®^ ^E^ 0, (3.37) 

where 0{—m) denotes the sheaf obtained by pulling back the corresponding 
sheaf over P-^ to T. In the case of bundles, the exact sequence (3.37) just means 
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that the dual space to the space of eigenvectors over the spectral curve is given 
by the cokernel of the linear map defined by >C^(A) — zir. A more complete 
discussion of the significance of this construction may be found in [AHH2], 
[AHH3]. 

3.4 Spectral Darboux Coordinates 

In this section we give a method for constructing the appropriate Darboux coor- 
dinates naturally associated to the spectral data discussed above, in which the 
Hamiltonians in the spectral ring IXks determine a Liouville generating func- 
tion in completely separated form. First, we shall give a purely computational 
description of these coordinates in terms of simultaneous solutions of polynomial 
equations. The significance of this prescription in terms of the eigenvector line 
bundles of the preceding section will follow. 

Let 

M(A,C):=^-ar, (3.38) 

and denote by M(A, Q the transpose of the matrix of cofactors. Then, over 
the spectral curve defined by the characteristic equation (3.34), the columns of 
M(A, are the eigenvectors of C{X) (or A/'(A)), and hence, generically, these are 
all proportional; i.e. M{\,Q has rank 1. Let Vq G be a fixed vector, and 
denote the solutions to the system of polynomial equations 

M(A, C)Vb = 0, Vb e C (3.39) 

as {A^, Cfj.}iJ.=i,...- Note that, due to the rank condition, there really are only two 
independent equations here, the other r — 2 following as linear consequences. 

The significance of these equations in relation to the spectral data is quite 
simple; they are the conditions that a section of the dual eigenvector line bundle 
£■ — > C should vanish. The solutions give the zeros of the components of the 
eigenvector determined by the vector Vq. As is well known in algebraic geometry, 
giving the divisor of zeros of any section of a line bundle amounts to giving the 
linear equivalence class of the bundle itself. It follows, since the bundle E ^ C 
is of degree g + r — 1, that there will in general be 5' -|- r — 1 zeros. However this 
is not necessarily the number of solutions to (3.39), since some of the zeros may 
be over A = 00. In fact, we may distinguish two cases of particular interest as 
follows. In order to characterize the spectrum over A = 00, define 



C{X) :=£(A)/A^. 



(3.40) 
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Case (i). Vq is an eigenvector of £(oo). In this case, r — 1 of the zeros are 
over A = oo (the only point omitted over oo being the one corresponding to the 
eigenvalue of Vq.) Hence, there are only g finite solutions pairs {A^, C^} ~, 
and these are generically independent, when viewed as functions on the phase 
space Ofj-^ (case (a)) or O^To (case (b)). Moreover, they are invariant under 
the action of the reduction group for both cases, since this leaves the space 
[Vq] invariant, and hence they project to functions on the reduced space O^j^^ 
(case (a)) or Ojj-^'^'^ (case (b)). In view of the dimensions given in the Table of 
Dimensions, Section 3.2, the projected functions provide coordinate systems on 
the reduced spaces. On the prereduced spaces, we must supplement these with 
a further r — 1 pairs of coordinate functions, which we define as follows. Choose 
a basis where Lq (case (a)) or Y (case (b)) is diagonal, and Vq — {1,0 . . . , 0)^. 
Then let 



ln(Li)ii + I YIj^2,3i^i ^"^'<P^ - Pj) case (a) 
ln(Lo)ii for case (b) (3.41a) 

Pi:=(Lo)ii, z = 2,...r. (3.41b) 

The pairs {g^, Pi}i=2,...r provide the remaining coordinates required. 

Case (ii). Vq is not an eigenvector of >C(oo) and, furthermore, Vq ^ Im (>C(oo) — 
ZiT) for any eigenvalue Si(oo) of £(oo). In this case, none of the zeros are 
over A = oo, and there are generically g + r — 1 independent solution pairs 
{^fj,iCij} ^+r-i (^-^Q)- These then provide a coordinate system on 

the prereduced space Ofj-^ (case (a)) or C^vb (case (b)). 

We then have the following fundamental result. 

Theorem 3.2. 1. I/Vq ^ Im (>C(oo) — 2^11) for any eigenvalue ii(oo) of C{oo), 
the functions {X^Xf^} I g+r-i define a Darboux coordinate system on Ofj-^ (case 
(a) or Oj^fg (case (b)). The orbital symplectic form is therefore: 

g+r-l 

J2 dXf^AdCf,. (3.42a) 

2. IfVo is an eigenvector of Lq (case (a) orY (case (b)), the functions {Xf^X^} i ^ 
project to Darboux coordinates on the reduced spaces Oj^^ (case (a)) or Ojj-^'^'^ 
(case (b)), so the reduced orbital symplectic form is: 



dX^ A c^Cm- (3-42b) 
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3. If Vq is an eigenvector of Lq (case (a)) or Y (case (h)), the functions 
{X^, Cfx: QiT Pi} 1^=1 i=2 r ^^fi"^^ ^ Darhoux coordinate system on Of^^ (case 
(a) or Cjv'o (case (h)), so the orbital symplectic form is: 

9 r 

t^orb ^^d\^^ dCi^ + ^ % A dPi. (3.42c) 

In the following section, we consider a number of elementary examples of 
the above theorem. We shall see that the resulting spectral Darhoux coordinates 
do, indeed, generalize the hyperellipsoidal coordinates that were encountered in 
the examples of Section 1. 

3.5 Examples 

We begin by considering the simplest possible case; namely, where Ao(A) has 
only one pole, at A = cti, and r = 2 or 3. This just corresponds to coadjoint 
orbits of the finite dimensional Lie algebras s[(2) and s[(3). Then we consider the 
case sl(2, R)"*" for arbitrary n, with rank(A^i) = ki = 1 for alH = 1, . . .n. This 
reproduces the hyperellipsoidal coordinates for the finite dimensional examples 
of Sees. 1.1 and 2.5.1 (cf. [Mo], such as the Neumann oscillator. Finally, we 
consider the case su(l, 1)"*", which provides the appropriate complex coordinates 
for the nonlinear Schrodinger equation, as discussed in Sees. 1.2 and 2.5.2. 

(a) Single poles: n = 1 

(a.l) Take g = s[(2,R), and (without loss of generality), ai = 0. Then the 
dimension of a generic orbit is dim — 2. We parametrize Ao(A) as follows: 

and choose 

y:-{l V„-{1). (3.44) 

The characteristic equation is then 

det {C{X) - zir) = z^ -a^ -ur = ^. (3.45) 

In this case, Vq is an eigenvector of Y and the genus of the spectral curve is = 0, 
so there are no {Ap,C/i}'s. The single pair of spectral Darboiix coordinates is 
thus 

= In u, P2 = a. (3.46) 
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It is easily verified that, relative to the Lie Poisson structure, they satisfy 

{92,P2} = 1. (3.47) 
(a.2) We consider the same orbit as in (a.l), but choose 

Y:=[1 I). (3.48) 

In this case, Vq is not an eigenvector of Y. The genus is still g = 0, but the 
equation (3.39) now has a finite solution, giving the Darboux coordinate pair 

Xi = -u, Ci = --- (3.49) 
u 

These are verified to also satisfy 

{Ai,Ci} = 0. (3.50) 

(a.3) Take g = s[(3, R), and again, cti = 0. Then the dimension of a generic 
orbit with n = 1 is dim C^vb = 6. We parametrize N'o{X) as: 



J\fo{X) = iVi := u a e , (3.51) 




and choose 




y= 1 , . (3.52) 




Again, Vq is an eigenvector of Y, but the spectral curve has genus g = 1, and is 
realized as a 3-fold branched cover of . We therefore find one Darboux coordi- 
nate pair (Ai,(^i), corresponding to a finite zero of the eigenvector components, 
plus two further pairs, (g2, P2, Qs, P3), corresponding to zeros over A = 00: 

. 1 /, ev uf\ ^ uva + uvb — ev^ — fu^ 
Ai = 7^M>-o + — , Ci = 



2 \ u V ) ^ —uva + uvh — ev^ + ju? 

q2 = In u, qs = In v, P2 = a, P3 = h. (3.53) 

Again, it is easily verified directly that these form a Darboux system, with 
nonvanishing Lie Poisson brackets 



{Ai,Ci} = l, {52,^2} = !, {53,^3} = !. 



(3.54) 
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(b) Now consider the case 0+ = s/(2,R)+, with arbitrary n, but ranks ki = 1 
for all i, and hence det(A^i) = for all the residue matrices A^^. For general Y, 
Af{X) then has the form 



^f{X) = A n +77 2' , (3.55) 

where {x^, yi}i=i,...n form a Darboux system on the reduced Moscr space, which 
is identified with M."^"^ / {^2)^ ■ In this case, the characteristic equation defining 
the invariant spectral curve C is 

det(>C(A) - zh) = -2^ + o(A)P(A) = 0, (3.56) 

where 

P(A) = -(a^ + &c)A" + Pn-iA"-' + . . . . (3.57) 

In particular, this gives eq. (1.34) for the case a = c = 0, b = — |. Thus, C is 
hyperelliptic, a 2-sheeted branched cover of P^, with 2n — 1 or 2n finite branch 
points, depending on whether or not + be vanishes. The genus is therefore 
generically ]j = n — 1. The dimension of the coadjoint orbit is dim Oj\/q = 2n. 
The matrix M{X, () is 

M= \ "^+2 2.i=l A-a, ^ 2.^=1 \ 



Taking 



Fo = ( J ) , (3.59) 



if c 7^ 0, Vo is not an eigenvector of Y , so the full set of n spectral Darboiix 
coordinate pairs {A^, CAt}At=i,...n are given by: 



n 9 



y - — — + 2c = 0, (3.60a) 



1 " 

Cm = -«+^Et^^' ^3.60b) 
fj, = 1, . . .n. 

These are therefore hyperelliptic coodinates {A^} and their conjugate momenta 
{Cfj,}- In the case c = 0, Vq is an eigenvector of Y, and eqs. (3.60a,b) are replaced 

by 
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yielding only n — 1 pairs of Darbonx coordinates {A^, Cfj.}iJ,=i,...n-i, since one of 
the zeros of the eigenvector components lies over A = oo. We must therefore 
complete the system by defining the additional pair 



^ n 1 

9:=ln(-J]x2), P■.^-J2x^y^. (3.62) 
It is easily verified directly that 



'2^ 2 

i=i i=i 



u = -de, (3.63) 



where 



\Tr-J^C,d\, + Pdq ifc=0. 

(c) NLS Equation: iu(l,l)+ 

Taking the orbit Oj^^ in iu(l, l)"*"* as parametrized in eqs. (2.54), (2.55), with 
y = 0, the symplectic submanifold Off^ C Oj^^ is defined by the constraint 

n 

E = 0- (3-65) 

The spectral Darboux coordinates {q, P, A^, Cm}m=1v Ti-i ^^"^ then given by eqs. 
(1.52a-c). It is easily verified in this case that the orbital symplectic form re- 
stricted to Ofj-^ is 

cjorb = —d9, (3.66) 

where 

n n—2 

so {q, P, A^, CiJ.}iJ.=i,...n-2 do, indeed, define a Darboiix coordinate system. 

A similar construction holds for the case of the sine-Gordon equation (Sec- 
tion 1.3), where the relevant algebra is the twisted loop algebra su(2)+, obtained 
by a suitable combination of discrete and continuous reductions. The orbits are 
parametrized by eqs. (1.60), (1.61a,b), and the relevant spectral Darboux coor- 
dinates determined by eqs. (1.84a,b). Details may be found in [HW]. 

In the last section we explain how, in the general case, these spectral 
Darboux coordinates lead directly to a linearization of the AKS fiows through 
the Liouville- Arnold integration procedure. In each case the relevant linearizing 
map turns out to be the Abel map to the Jacobi variety of the spectral curve. 
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3.6 Liouville- Arnold Integration 

Using the spectral Darboux coordianates, we may define the local equivalent of 
the "canonical" 1-form 

g r 

e:=Y^ C^^dX^ + PidQi (3.68a) 

/u=l i=2 
g r 

where 

•= «(-^a*)Cm- (3-69) 

(Note that for the examples given above, this actually is the canonical 1-form 
on M2n/(^2)n^ C2"/(Z2)", vlcwcd as the cotangent bundle of W/{'L2Y and 
C"^/ (Z2)", respectively.) On the Liouville- Arnold torus T, defined by taking the 
simultaneous level sets 

Pia — Cia, Pi — Ci (3.70) 
of the spectral invariants, we have 

6'It = dS{Xi,...X~,q2,...qr,Pia,Pi), (3.71) 

where S'(Ai, . . . , A~, g2, • • • , Qr, Pia, Pi) is the Liouville generating function to the 
canonical coordinates conjugate to the invariants {Pia, Pi)- Eq. (3.71) can be in- 
tegrated by viewing z = z{X, Pia, Pi) as a, meromorphic function on the Riemann 
surface of the spectral curve C. 

SiX^, Qi, Pra, Pi) = r '-^^"^^dX + Y ^iPi^ (3-72) 

where = z^^{X^, Pia, Pi) is essentially determined implicitly by the spectral 
equation 

V{X^,z^{X^,Pia,Pi)) = Q. (3.73) 

The linearizing coordinates for AKS flows are then given, as usual, by differen- 
tiation of S with respect to the invariants: 

dS 4^ f^^ 1 dz „ dh ^ , 



^ dS f^^ 1 dz „ dh 

= 'aB' = 1^ / T^'^dX + qi = -^t + Qi^o, 

i^JK (3.74b) 
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where, from the exphcit structure (3.27) of the characteristic polynomial given 
in Proposition 3.1, we obtain, by implicit differentiation, that the integrands of 
eqs. (3.74a,b) are of the form 



where 



(-Pi Pi) Z]2<ii<i2-"<ij_27^i Pii ■ • ■ Pij-2 

and e = for case (a) 
Rij:= { ^ ^ (3.76) 

and 6=1 for case (b) . 

The point to note is that the differentials {uia}, {uJi} appearing in eqs. (3.75a,b) 
are, respectively, abelian differentials of the first and third kinds on the Riemann 
surface defined by C, the latter having their poles at the points (ooi, . . . ,oOj.) 
over A = oo. 



Theorem 3.3 [AHH3]. Theg differentials {uJia}i=i,...g in eq. (3.75a) form a 
basis for the space H'^{S, K~) of abelian differentials of the first kind (where K~ 
denotes the canonical bundle). The linear flow equation (3.74a) may therefore 
be expressed as: 

A(V) = B + Ut, (3.77) 

where A : S^C — > C^/F is the Abel map, and B, U G are obtained by 
applying the inverse of the g xg normalizing matrix M, with elements 

M^,{ia) / ^ia, (3.78) 

to the vectors C,H e with components Cia and respectively (the pair 

(ia) viewed as a single coordinate label in ). The r—1 differentials {u!i}i=2,...r 
in eq. (3.75b) are abelian differentials of the third kind with simple poles at ooj 
and ooi, and residues +1 and —1, respectively. 

Comparing this general formula with the specific cases (1.23), (1.55), (1.56), 
(1.88a,b), for the examples of Section 1, we see that this provides the generaliza- 
tion that was required, expressing all linearized AKS fiows on rational coadjoint 
orbits of sl{r)+* and its reductions through the Abel map. 
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It is possible, moreover, to invert the map expressing these flows, by ex- 
pressing any symmetric function of the coordinates {A^, Cm}^=i ^ terms of 
the Riemann theta function associated to the curve C. For example, in view of 
eq. (3.74b), the coordinates {qi}i=2,...r themselves are expressed as such symmet- 
ric functions through abelian integrals of the third kind. Applying the reciprocity 
theorem relating the two kinds of abelian integrals {viz. [AHH3]), we obtain: 



Corollary 3.4 [AHH3]. For a suitable choice of constants {ci, fi}i=2, 
coordinate functions {qi{t)} satisfying eq.(3.74b) are given by: 



the 



qi{t) = In 



6>(B + tU - A(ooj) - K) 
^(B + tU- A(ooi) -K) 



+ eit + fi, 



(3.79) 



where K e is the Riemann constant. 



This generalizes the theta function formula (1.57) giving the solution of the 
NLS equation. Similar formulae exist e.g., for the sine-Gordon equation ([HW]) 
and many other systems that can be cast in terms of commuting AKS flows in 
rational coadjoint orbits of loop algebras. Aside from technical complications 
resulting, e.g., from the reduction procedure or the imposition of further sym- 
plectic constraints, or from the presence of further singularities in the spectral 
curve, the procedure is largely algorithmic. It provides a very general setting 
for the explicit application of the Liouville-Arnold integration procedure to a 
wide class of known - and yet to be discovered - integrable Hamiltonian systems. 
Moreover, the moment map embedding method makes it possible to treat both 
the intrinsically flnite dimensional systems, and those systems corresponding to 
flnite dimensional sectors of integrable systems of PDE's (solitons, flnite band 
solutions, etc.) on exactly the same footing. 
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